Quantitative test of the barrier nucleosome model for statistical positioning of nucleosomes up- and downstream of transcription start sites

The positions of nucleosomes in eukaryotic genomes determine which parts of the DNA sequence are readily accessible for regulatory proteins and which are not. Genome-wide maps of nucleosome positions have revealed a salient pattern around transcription start sites, involving a nucleosome-free region (NFR) flanked by a pronounced periodic pattern in the average nucleosome density. While the periodic pattern clearly reflects well-positioned nucleosomes, the positioning mechanism is less clear. A recent experimental study by Mavrich et al. argued that the pattern observed in S. cerevisiae is qualitatively consistent with a `barrier nucleosome model', in which the oscillatory pattern is created by the statistical positioning mechanism of Kornberg and Stryer. On the other hand, there is clear evidence for intrinsic sequence preferences of nucleosomes, and it is unclear to what extent these sequence preferences affect the observed pattern. To test the barrier nucleosome model, we quantitatively analyze yeast nucleosome positioning data both up- and downstream from NFRs. Our analysis is based on the Tonks model of statistical physics which quantifies the interplay between the excluded-volume interaction of nucleosomes and their positional entropy. We find that although the typical patterns on the two sides of the NFR are different, they are both quantitatively described by the same physical model, with the same parameters, but different boundary conditions. The inferred boundary conditions suggest that the first nucleosome downstream from the NFR (the +1 nucleosome) is typically directly positioned while the first nucleosome upstream is statistically positioned via a nucleosome-repelling DNA region. These boundary conditions, which can be locally encoded into the genome sequence, significantly shape the statistical distribution of nucleosomes over a range of up to ~1000 bp to each side.Comment: includes supporting materia

Inheritance of protection from osmotic stress

Exposure of mother worms to mild osmotic stress induces gene expression changes in offspring that protect them from strong osmotic stress. Inheritance of protection is now shown to depend on altered insulin-like signalling in the maternal germline, which confers protection through increased expression of zygotic gpdh-2, a rate-limiting enzyme in glycerol biosynthesis

Replication and Active Demethylation Represent Partially Overlapping Mechanisms for Erasure of H3K4me3 in Budding Yeast

Histone modifications affect DNA–templated processes ranging from transcription to genomic replication. In this study, we examine the cell cycle dynamics of the trimethylated form of histone H3 lysine 4 (H3K4me3), a mark of active chromatin that is viewed as “long-lived” and that is involved in memory during cell state inheritance in metazoans. We synchronized yeast using two different protocols, then followed H3K4me3 patterns as yeast passed through subsequent cell cycles. While most H3K4me3 patterns were conserved from one generation to the next, we found that methylation patterns induced by alpha factor or high temperature were erased within one cell cycle, during S phase. Early-replicating regions were erased before late-replicating regions, implicating replication in H3K4me3 loss. However, nearly complete H3K4me3 erasure occurred at the majority of loci even when replication was prevented, suggesting that most erasure results from an active process. Indeed, deletion of the demethylase Jhd2 slowed erasure at most loci. Together, these results indicate overlapping roles for passive dilution and active enzymatic demethylation in erasing ancestral histone methylation states in yeast

Quantitative Analysis of Histone Modifications: Formaldehyde Is a Source of Pathological N6-Formyllysine That Is Refractory to Histone Deacetylases

Aberrant protein modifications play an important role in the pathophysiology of many human diseases, in terms of both dysfunction of physiological modifications and the formation of pathological modifications by reaction of proteins with endogenous electrophiles. Recent studies have identified a chemical homolog of lysine acetylation, N[superscript 6]-formyllysine, as an abundant modification of histone and chromatin proteins, one possible source of which is the reaction of lysine with 3′-formylphosphate residues from DNA oxidation. Using a new liquid chromatography-coupled to tandem mass spectrometry method to quantify all N[superscript 6]-methyl-, -acetyl- and -formyl-lysine modifications, we now report that endogenous formaldehyde is a major source of N[superscript 6]-formyllysine and that this adduct is widespread among cellular proteins in all compartments. N[superscript 6]-formyllysine was evenly distributed among different classes of histone proteins from human TK6 cells at 1–4 modifications per 10[superscript 4] lysines, which contrasted strongly with lysine acetylation and mono-, di-, and tri-methylation levels of 1.5-380, 5-870, 0-1400, and 0-390 per 10[superscript 4] lysines, respectively. While isotope labeling studies revealed that lysine demethylation is not a source of N[superscript 6]-formyllysine in histones, formaldehyde exposure was observed to cause a dose-dependent increase in N[superscript 6]-formyllysine, with use of [[superscript 13]C,[superscript 2]H[subscript 2]]-formaldehyde revealing unchanged levels of adducts derived from endogenous sources. Inhibitors of class I and class II histone deacetylases did not affect the levels of N[superscript 6]-formyllysine in TK6 cells, and the class III histone deacetylase, SIRT1, had minimal activity (<10%) with a peptide substrate containing the formyl adduct. These data suggest that N[superscript 6]-formyllysine is refractory to removal by histone deacetylases, which supports the idea that this abundant protein modification could interfere with normal regulation of gene expression if it arises at conserved sites of physiological protein secondary modification

Discovery and Annotation of Functional Chromatin Signatures in the Human Genome

Transcriptional regulation in human cells is a complex process involving a multitude of regulatory elements encoded by the genome. Recent studies have shown that distinct chromatin signatures mark a variety of functional genomic elements and that subtle variations of these signatures mark elements with different functions. To identify novel chromatin signatures in the human genome, we apply a de novo pattern-finding algorithm to genome-wide maps of histone modifications. We recover previously known chromatin signatures associated with promoters and enhancers. We also observe several chromatin signatures with strong enrichment of H3K36me3 marking exons. Closer examination reveals that H3K36me3 is found on well-positioned nucleosomes at exon 5′ ends, and that this modification is a global mark of exon expression that also correlates with alternative splicing. Additionally, we observe strong enrichment of H2BK5me1 and H4K20me1 at highly expressed exons near the 5′ end, in contrast to the opposite distribution of H3K36me3-marked exons. Finally, we also recover frequently occurring chromatin signatures displaying enrichment of repressive histone modifications. These signatures mark distinct repeat sequences and are associated with distinct modes of gene repression. Together, these results highlight the rich information embedded in the human epigenome and underscore its value in studying gene regulation

Inheritance of Acquired Behaviour Adaptations and Brain Gene Expression in Chickens

Background: Environmental challenges may affect both the exposed individuals and their offspring. We investigated possible adaptive aspects of such cross-generation transmissions, and hypothesized that chronic unpredictable food access would cause chickens to show a more conservative feeding strategy and to be more dominant, and that these adaptations would be transmitted to the offspring. Methodology/Principal Findings: Parents were raised in an unpredictable (UL) or in predictable diurnal light rhythm (PL, 12:12 h light:dark). In a foraging test, UL birds pecked more at freely available, rather than at hidden and more attractive food, compared to birds from the PL group. Female offspring of UL birds, raised in predictable light conditions without parental contact, showed a similar foraging behavior, differing from offspring of PL birds. Furthermore, adult offspring of UL birds performed more food pecks in a dominance test, showed a higher preference for high energy food, survived better, and were heavier than offspring of PL parents. Using cDNA microarrays, we found that the differential brain gene expression caused by the challenge was mirrored in the offspring. In particular, several immunoglobulin genes seemed to be affected similarly in both UL parents and their offspring. Estradiol levels were significantly higher in egg yolk from UL birds, suggesting one possible mechanism for these effects. Conclusions/Significance: Our findings suggest that unpredictable food access caused seemingly adaptive responses in feeding behavior, which may have been transmitted to the offspring by means of epigenetic mechanisms, including regulation of immune genes. This may have prepared the offspring for coping with an unpredictable environment. Citation: Nätt D, Lindqvist N, Stranneheim H, Lundeberg J, Torjesen PA, et al. (2009) Inheritance of Acquired Behaviour Adaptations and Brain Gene Expression in Chickens. PLoS ONE 4(7): e6405. doi:10.1371/journal.pone.0006405 Editor: Tom Pizzari, University of Oxford, United Kingdom Received: March 26, 2009; Accepted: June 30, 2009; Published: July 28, 2009 Copyright: © 2009 Nätt et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: This project was funded by the Swedish Research Council (VR; www.vr.se; grant nrs 50280101 and 50280102) and the Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning (Formas; www.formas.se; grant no 221-2005-270). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the mauscript. Competing interests: The authors have declared that no competing interests exist.  Original Publication:Daniel Nätt, Niclas Lindqvist, Henrik Stranneheim, Joakim Lundeberg, Peter A. Torjesen and Per Jensen, Inheritance of Acquired Behaviour Adaptions and Brain Gene Expression in Chickens, 2009, PLoS ONE, (4), 7, e6405.http://dx.doi.org/10.1371/journal.pone.0006405Copyright: Author

Exploring the Contextual Sensitivity of Factors that Determine Cell-to-Cell Variability in Receptor-Mediated Apoptosis

Stochastic fluctuations in gene expression give rise to cell-to-cell variability in protein levels which can potentially cause variability in cellular phenotype. For TRAIL (TNF-related apoptosis-inducing ligand) variability manifests itself as dramatic differences in the time between ligand exposure and the sudden activation of the effector caspases that kill cells. However, the contribution of individual proteins to phenotypic variability has not been explored in detail. In this paper we use feature-based sensitivity analysis as a means to estimate the impact of variation in key apoptosis regulators on variability in the dynamics of cell death. We use Monte Carlo sampling from measured protein concentration distributions in combination with a previously validated ordinary differential equation model of apoptosis to simulate the dynamics of receptor-mediated apoptosis. We find that variation in the concentrations of some proteins matters much more than variation in others and that precisely which proteins matter depends both on the concentrations of other proteins and on whether correlations in protein levels are taken into account. A prediction from simulation that we confirm experimentally is that variability in fate is sensitive to even small increases in the levels of Bcl-2. We also show that sensitivity to Bcl-2 levels is itself sensitive to the levels of interacting proteins. The contextual dependency is implicit in the mathematical formulation of sensitivity, but our data show that it is also important for biologically relevant parameter values. Our work provides a conceptual and practical means to study and understand the impact of cell-to-cell variability in protein expression levels on cell fate using deterministic models and sampling from parameter distributions

A complete set of nascent transcription rates for yeast genes

The amount of mRNA in a cell is the result of two opposite reactions: transcription and mRNA degradation. These reactions are governed by kinetics laws, and the most regulated step for many genes is the transcription rate. The transcription rate, which is assumed to be exercised mainly at the RNA polymerase recruitment level, can be calculated using the RNA polymerase densities determined either by run-on or immunoprecipitation using specific antibodies. The yeast Saccharomyces cerevisiae is the ideal model organism to generate a complete set of nascent transcription rates that will prove useful for many gene regulation studies. By combining genomic data from both the GRO (Genomic Run-on) and the RNA pol ChIP-on-chip methods we generated a new, more accurate nascent transcription rate dataset. By comparing this dataset with the indirect ones obtained from the mRNA stabilities and mRNA amount datasets, we are able to obtain biological information about posttranscriptional regulation processes and a genomic snapshot of the location of the active transcriptional machinery. We have obtained nascent transcription rates for 4,670 yeast genes. The median RNA polymerase II density in the genes is 0.078 molecules/kb, which corresponds to an average of 0.096 molecules/gene. Most genes have transcription rates of between 2 and 30 mRNAs/hour and less than 1% of yeast genes have >1 RNA polymerase molecule/gene. Histone and ribosomal protein genes are the highest transcribed groups of genes and other than these exceptions the transcription of genes is an infrequent phenomenon in a yeast cell

Small RNA-Directed Epigenetic Natural Variation in Arabidopsis thaliana

Progress in epigenetics has revealed mechanisms that can heritably regulate gene function independent of genetic alterations. Nevertheless, little is known about the role of epigenetics in evolution. This is due in part to scant data on epigenetic variation among natural populations. In plants, small interfering RNA (siRNA) is involved in both the initiation and maintenance of gene silencing by directing DNA methylation and/or histone methylation. Here, we report that, in the model plant Arabidopsis thaliana, a cluster of ∼24 nt siRNAs found at high levels in the ecotype Landsberg erecta (Ler) could direct DNA methylation and heterochromatinization at a hAT element adjacent to the promoter of FLOWERING LOCUS C (FLC), a major repressor of flowering, whereas the same hAT element in ecotype Columbia (Col) with almost identical DNA sequence, generates a set of low abundance siRNAs that do not direct these activities. We have called this hAT element MPF for Methylated region near Promoter of FLC, although de novo methylation triggered by an inverted repeat transgene at this region in Col does not alter its FLC expression. DNA methylation of the Ler allele MPF is dependent on genes in known silencing pathways, and such methylation is transmissible to Col by genetic crosses, although with varying degrees of penetrance. A genome-wide comparison of Ler and Col small RNAs identified at least 68 loci matched by a significant level of ∼24 nt siRNAs present specifically in Ler but not Col, where nearly half of the loci are related to repeat or TE sequences. Methylation analysis revealed that 88% of the examined loci (37 out of 42) were specifically methylated in Ler but not Col, suggesting that small RNA can direct epigenetic differences between two closely related Arabidopsis ecotypes

Missing value imputation for microarray gene expression data using histone acetylation information

<p>Abstract</p> <p>Background</p> <p>It is an important pre-processing step to accurately estimate missing values in microarray data, because complete datasets are required in numerous expression profile analysis in bioinformatics. Although several methods have been suggested, their performances are not satisfactory for datasets with high missing percentages.</p> <p>Results</p> <p>The paper explores the feasibility of doing missing value imputation with the help of gene regulatory mechanism. An imputation framework called histone acetylation information aided imputation method (HAIimpute method) is presented. It incorporates the histone acetylation information into the conventional KNN(<it>k</it>-nearest neighbor) and LLS(local least square) imputation algorithms for final prediction of the missing values. The experimental results indicated that the use of acetylation information can provide significant improvements in microarray imputation accuracy. The HAIimpute methods consistently improve the widely used methods such as KNN and LLS in terms of normalized root mean squared error (NRMSE). Meanwhile, the genes imputed by HAIimpute methods are more correlated with the original complete genes in terms of Pearson correlation coefficients. Furthermore, the proposed methods also outperform GOimpute, which is one of the existing related methods that use the functional similarity as the external information.</p> <p>Conclusion</p> <p>We demonstrated that the using of histone acetylation information could greatly improve the performance of the imputation especially at high missing percentages. This idea can be generalized to various imputation methods to facilitate the performance. Moreover, with more knowledge accumulated on gene regulatory mechanism in addition to histone acetylation, the performance of our approach can be further improved and verified.</p